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FOREWORD 


This report wes prepared by Captain George P. Frost and 
Dr. Arthur J. Wenneistrom of the Fluid Dynemics Facilities 
Research Laboratory, Aerospace Research Laboratories, Wright-~ 
Patterson Air Force Base, Ohio. 


The report presents results from a portion of the effort of 
the Fluid Nachinery Research Group supervised by Dr. Arthur J. 
Wennerstron and was conducted under Work Unit 09 of Project 7665, 
"Aerospace Simulation Techniques Research" under the overall 
direction of Mr. Elmer G. Johnson. 


oe. 
oe 


RN ae op lll ede <A Il Mee AI HY ra, aie acim csemeteedinnes hl mys 


A 


a am wet ow 


ne a i eel el Ree lente anes ee hh i el 


ee 
“ 


ee 


——er 


ABSTRACT 


This report describes a technique which has been developed 
for use in the design of axial compressor airfoils with camber 
lines of arbitrary shape. The slope of the camber line at 
several points on a streamsurface is determined from the air 
angles at thesé points as well as the incidence and deviation 
angle distributions for the blade. A camber line is produced by 
Fitting a smooth curve segment through each pair of points from 
the leading to the trailing edge. A thickness distribution is 
applied to this camber line to produce the blade element. A 
computer program which uses this technique to produce blade 
elements. stack them, and then determine coordinates for plane 
surfaces through the resultant blade is aiso described. 
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SECTION I 
INTRODUCTION 


The traditional approach te the selection of blades for 
axial-flow compressors has been to choose a specific type of 
airfoil, subsequently to adjust the parameters defining that 
airfoil to suit prevailing aerodynamic conditions. Many 
successful designs have been accomplished in this manrer using 
such blades as NACA 65-series, double circviar arc, multiple 
circular arc, etc. This approach has the virtue of employing 
airfoil sections about which something is usually known with 
respect to losses, deviation angle, and operating limits. 
Reference 1, an earlier publication of this laboratory, is 
representative of this approach. 


A typical contemporary design will be accomplished using a 
streamline curvature or matrix through-flow analysis incorporat- 
ing computing stations internal to blade rows as well as at blade 
edges and in free spaces. Some criterion for optimization of 
Dlade shape must be chosen. (The authors of this report choose 
to achieve a particular shape of static pressure distribution 
along streamlines.) The spanwise distributions of inlet and 
outlet relative flow angles, loss coefficients, etc. will have 
been determined through some preliminary design procedure. A 
rule for determining deviation angle and distributing it along 
streamlines will have been selected. Subsequently, the aero- 
dynamic blade design procedure consists of assuming a blade 
geometry, performing an aerodynamic analysis using specified 
relative flow angles as input data, and then repeating this 
procedure as many times as necessary, varying the parameters 
describing the airfoil sections, until the result of the aero- 
dynamic analysis adequately satisfies the chosen optimization 
criterion. 


There are three principal disadvantages to employing this 
procedure for the design of transonic and supersonic blade rows. 
First, the iteration required is time consuming and, to some 
degree, laborious. Many adjustments to the geometry specified 
may be required before optimization objectives are met. 

Although the bulk of the calculations are performed by computer, 
the designer generally must still examine and evaluate the 
aerodynamic result and decide what to change, and how much, for 
the next attempt. Second, when airfoil shape is restricted to 
any particular class of airfoil, it will rarely be possible to 
achieve the optimization objectives as closely as might be 
desired, everywhere along the span. Some compromise will nearly 
always be necessary. Third, streamline curvature and matrix 
through-flow analyses sometimes experience difficulty in finding 
solutions for certain combinations of high relative Mach number 
and high absolute Mach number near choking when relative flow 
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angle is specified as input data within blade rows. This is a 
numerical problem related to the construction of these computer 
programs and can vary considerably from one program to another. 


A design method which eliminates or appreciably reduces 
these problems consists of specifying total temperature (in 
rotors) or the product of radius and whirl velocity (in stators) 
as input data to the aerodynamic analysis program and then 
fitting airfoils of arbitrary shape to the resulting relative 
flow angies. In so doing, one loses the data base which might 
be associated with a particular class cf airfoil, but the value 
of such a geometrically-related data base is questionable for 
transonic and supersonic sections. Foliowing this approach, 
the aerodynamic analysis can be optimized with considerably 
fewer iterations than are usually required with specified 
geometry. In most cases, optimization objectives can be achieved 
on aearly every streamsurface. Some interaction is required 
between the blade design program and the aerodynamic analysis 
to insure that the blade lean angles and blockages are kept 
up-to-date. However, the over-all effort required to achieve 
a satisfactory design using this method has been found to be 
substantially less. 


The principal difficulty in developing a procedure to 
define arbitrary airfoils consists of arriving at a method 
which produces aerodynamically attractive shapes which are in 
addition practical from a structural and manufacturing viewpoint. 
This report describes one such method developed at the Aerospace 
Research Laboratories. The method has been incorporated into 
a computer program which is an extension of the work reported 
in Reference 1. In additio., to determining the shape of the 
airfoils on aerodynamic surfaces, section properties are 
computed, the blade is stacked, and Cartesian coordinates are 
determined for manufacturing purposes. 


The overali design technique is described in Section II. 
The mathematical details of implementing the technique in a 
calculation procedure are described in Sections IIT and IV. A 
method of producing the optimal camber line on a streamsurface 
is treated in Section III; the calculations related to other 
aspects of the blade design procedure are discussed briefly in 
Section IV. Sections V, VI, and VII present the details of and 
use of a computer program which incorporates this design 
technique, currently used at the Aerospace Research Laboratories 
in the design vf axial compressor airfoils. 
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SECTION II 
DESCRIPTION OF THE TECHNIQUE 


1. TECHNIQUE OVERVIEW 


The technique described in this report uses an iterative 
procedure to produce an “optimal” camber line on each stream- 
surface. A thickness destribution is applied to each camber 
line, and the resulting blade elements are stacked tc produce 
the desired airfoil. The iechnique requires the designer to 
specify the incidence distribution radially at the blade Leading 
edge, the location of the stack axis and the stacking offsets 
of zach streamsurface section centroid therefrom, the pa‘ameters 
of the thickness distribution, and the chordwise distribution 
of deviation angle along the span. 


The optimization criterion which has been selected for the 
section camber line is to maximize the absolute value of the 
minimum radius of curvature on the camber line. The "optimal" 
camber line is chosen from a set of camber lines containing the 
minimum number of inflection points. This original set of 
camber lines is generated by varying the second derivative at 
the leading edge. 


The starting point of the design technigue is output data 
from an aerodynamic analysis of a particular blade row which 
has been generated by specifying a parameter other than blade 
geometry, such as those suggested in the Introduction, across 
the blade row. This data is in the form of the meridional 
coordinates of the streansurfaces and the chordwise distribution ) 
of the relative air angles on each streamsurface. The essential 
steps of the technique itself are described in a qualitative 
sense in the remainder of this section. 


2. DETAILS OF THE TECHNIQUE 


The first step of the overall procedure is to determine the 
optimal blade section on each streamsurface. This in itself 
is a multiple-step process which incorporates an iteration with 
solidity to establish a camber line for each of a range of values | 
of the second derivative at the leading edge, a search procedure 
to choose the "optimal" camber line, and the application of a 
thickness distribution to this camber line. 


The procedure begins with an initial estimate of solidity 
on the streamsurface. This estimate is made by applying the 
Stagger angle, assumed equal to the average of the inlet and i 
outlet relative air angles, to the meridional chord length, ! 
obtained by integrating along each streamline from the assumed 
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leading to trailing edge, to get a first estimate of the true 
chord. The solidity is then computed from this estimate, the 
mean streamsurface radius, aid the number of biades in the blade 
row. : 


The total deviation angle is computed from this estimate 
in some fashion, such as the modified Carter's Rule with 
appropriate constants. The deviation angle at each internal 
point is next determined as a designer-specified fraction of 
the total deviation. The required secticn angle at each 
internal point and at the trailing edge is then the difference 
‘between the relative air angle and the deviation angie, while 
the section angie at the leading edge is the difference between 


the relative air angle and the incidence angle. 


Each section camber Line is determined by fitting a segment . 
of a smooth curve, such as a cubic, between each pair of points 
from the leading to the trailing edge of the section. The slope 
at the endpoint of each segment matches the specified section 
angle there. The true chord length and the associated value of 
solidity can then be determined. If the solidity differs by 
more than a prescribed tolerance from the previous estimate, 
the steps subsequent to and including the total deviation 
determination are repeated for-the revised values of solidity 
until the desired tolerance is achieved. 


fhe entire procedure described thus far is repeated for a 
range of values of the second derivative at the leading edge. 
The resultant set of camber lines are inspected first to focus 
ohly on those which contain the minimum number of inflection 
points, and from these to choose the camber line which is 
“optimal” in the sense previously described. 


A thickness distribution is then applied to this camber 
line, and the procedure is repeated for each streamsurface. 
The blade sections are stacked, and the Cartesian coordinates 
of the resulting blade determined. In addition, the blade 
blockage and lean angle are computed at each appropriate strean- 
suxrface-computing station intersection point. 


As a final step, the designer inspects the resulting 
blockages, Cartesian centroid osfsets, blade lean angles, and 
coincidence cf the blade edges with the designated computing 
stations. If necessary, appropriate changes are made in the 
inputs to the aerodynamic analysis, and the entire procedure 
recycled until adequate overall coin: idence between the blade 
design and the aerodynamic analysis is achieved. 


SECTION Til 
THE SECTION CAMBER LIKE 


The First item required in the application of the technique 
described in the preceding section is the meridional chord length 
of the biade element, obtained by integrating along the strean- 
iine between the assumed ieading and trailing edges. This is 
accomplished by passing a spline curve through the meridional 
coordinates (x, vr) of the streamsurface. The slope of the 
streamsurface is calculated (as the slope of the spline-curve) 
at 100 points distributed uniformly on the x-axis (axially) 
between the edges of the biade section. The chord length, Cp, 
is obtained from the equation 


160 


(1) 


n=2 


An estimate of true chord is obtained by applying the 
stagger angle, assumed equal to the average of the inlet and 


outlet relative air angles, to the meridional chord so that the 
true chord estimate, C,, is 


BeetBre (2) 


The first estimate of solidity may be computed from the 
equation 


NC, 
* % [tastes] (3) 
i 9 


where N is the number of blades and ree, req are the radii of 
the streamsurface at the leading and trailing edges, 
respectively. 
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where 


The calculation of the deviation angle, 5, follows the 
NASA method (Reference 2, Equations 269 and 271) with an 
additional term, yY. 


_ 


_ n : a 
es: Soe ae i “Bb (8,.~i-6,.+8) saa (4) 


19 


is the variaticn from the reference deviation 

for a 10 percent-thick NACA 65-series thickness 
distribution 

is a correction factor for a blade shape witha a 
thickness distribution different from 2 65-serwies 
biade 


is a correction factor for a maximum thickness 
other than 10 percent 


is the slope of the deviation angle variation 
from reference deviation with camber 


is the solidity exponent (variablé with air 
inlet angle) 


is the relative air ang”™e at the particular edge 
is the incidence angle 


is the arbitrary extra deviation 


Solving Eq (4) for 6 yields 


For use 


specified by 


(5) 


in the calculation procedure, Kg_, 1, and Y are 
the designer. Several of the other quantities are 
ebtained from known quantities and figures of Reference 2: 

from Figure 161; K§+ » Figure 172; m, Figure 1663; and b, 
Figure 164. 
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The blade angle, a, is established at several points 
across the biade element from the relative air angle, 8, 
modified by the proper consideration of incidence or deviation. 
At the leading edge, 


ge = Bee - 1 (6) 


and at the trailing edge, 


Cte = Bie =e (7) 


At the other points, the biade angles are determined by 
subtracting a fracticn f of the traitling-edge deviation fron 
the relative air angle. At esch internal point j, 


57 F536 (8) 


The fraction fs is determined by radial interpoiation from the 
deviation distributians specified by the designer. 


The camber line if constructed by fitting a third order 
poiynemial (cubic) through each pair of points from the leading 
to the trailing edge of the section. Thus, each segment of the 
camper line is defined by equations of the form 


y = ax3 + bx? + cx + d (9) 
y' = 3ax2 + 2bx + ¢ | (30) 
y" = 6ax + 2b (11) 


As a result, each segment has at mest one inflection point 
{y" = 0) in its usefus range, and in most instances has none. 


The constants a, b, ¢c, and d for any particular segment can 
be expressed in terms of the endpoints (denoted by subscripts 
1 and 2) of that segment as 


(yo'-y¥:') - yi"(x2-%1) 
a= 19 
3 | (x22-x1?) = 2x1 (x2-x1) | are 


¥i- - 6aXx; 


= (23) 
c = y,' - 3axy? - 2bx, (14) 
d = yy - ax,? - bxy* - CX} (15) 
For simplicity, the leading edge of the camber line is 
pl.ced at the crigin of the coordinate system, resulting in 
the following boundary conditions for the first segment: 
At x=Q, y = 0 
y*' = tan dpp (16) 
) it 
yo = De 
At x = Xz, y’ = tan a, (17) 


With these boundary conditions, the apprcpriate values of 
a, b, c, and d for this segment are completely determined. 


The boundary conditions for the second and subsequent 
segments are specified at one endpoint by equating the first 
and second derivatives to the values for the preceding segment 
at the point of juncture; for example, for the second segment, 


= x = 
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(18) 
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and at the other endpoint, 
At x = X95 y’ = tan ao (19) 


From these conditions, a distinct set oF constants 
(a, b, ¢, d) are computed for the second segment. This same 
procedure is applied to each pair of points to preduce a 
camber line with continuous first and second derivatives all 
along its length. 


Mote that the first segment of the camber line reauires 
the specification of the second derivative y," at the leading 
edge. This boundary condition effects the constants of the 
first segment and thus the nature of this entire segment, includ- 
ing the conditions (y,, y,") at the other endpoint. Since the 
constants for the second segment are established From these 
conditions, and so on for the rest of the segments, the nature 
of the entire camber line depends on the value of y," specified 
at the leading edge. 


It has been found convenient to specify y," in terms of a 
non-dimensional parameter S/R,, the ratio of blade spacing, S, 
to the radius of curvature, Ry, at the leading edge. R, 2s given 


by the equation 
3 
[i + tan? tye |” 
R = (20) 


% i 
Yo 


and S$ is obtained fron 


ae Ld (21) 


5 ant, . Vas 
R. N | + tan2 age | (22) 
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Solving for y," gives 


aft + tan? a, |? « 
Te JOS (23) 
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which indicates that y," is the parameter S/Rp multipl 
constant. 


For @ particular vaiue of S/Kp, the true chord length of 
the resulting camber lines may be determined fron the endpoint 
of the final segment as 


C= Xe)? + (y,.)? (24) 


This value is used to compute a revised value of solidity, which 
is compared to the original estimate. If satisfactory 
coincidence has not been achieved, a corrected deviation angle 
is computed from the revised solidity, and the camber line 
reconstructed. This iteration is repeated until adequate 
coincidence has been obtained. 


Tt is difficult if not impossible to have an intvitive 
notion of a "good" value of yo" (hence, S/Ro). For some ranges 
of S/R 5, each camber line segment may contain en inflection 
point, while for other ranges, few if any segments may contain 
such a point. The authors of this report have assumed that, 
for aerodynamic as well as mechanical reasons, the most 
desirable airfoil among several matching the same flow angle 
at each computing station will Le one having the minimum number 
of inflection points between leading and trailing edge. This 
is accomplished by calculating the camber line for a broad range 
of values of S/R, and isolating the range in which the minimum 
total number of inflection points occurs. This range is examined 
with finer S/R, increments to generate a new set of camber lines 
on which the minimum absolute radius of curvature is identified. 
The value of 3/R, which produces the largest such radius is 
then made the mid-value of S/K, for the final search pass, 
using still finer S/R, increments. The camber line which 
possesses the largest value of the minimum radius of curvature 
at the conclusion of this search procedure is chosen as opt? 
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SECTION 1% 
OTHER ASPECTS OF THE CALCULATION PROCEDURE 
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Various other aspects of the calculation procedure are 
discussed briefly in this section. For greater detail on these 
topics, the reader's attention is directed to Reference 1, where 
these items are treated at some length as elements of the 
caiculation procedure ux which the subject procedure is a 
modification. 


1. THE SECTION THICKKESS DISTRIBUTION 


The thickness distribution which is applied to the camber 
line herein described is that referred to as the "Standard 
Thickness Distribution" in Reference 1. The distributicen is 
defined by two third-order polynomials, one from the leading edge 
to the point of maximum thickness, and another from there to tne 
trailing edge. At the point of juncture, the thickness end the 
first and second derivatives of thickness are equated. In order 
to prevent a reflex curvature from occurring in the thickness 
distribution near the leading edge, the second derivative of 
thickness is set equal to zero at the leading edge. The thick- 
ness of the leading and trailing edges is independently specified 
so that it need not be the same at both edges. At the leading 
edge, the blade surface is completed with a circular arc. At 
the trailing edge, the blade surface is truncated by connecting 
the two endpoints with a straight line. 


2. CARTESIAN COORDINATES FOR THE BLADE 


The preceding material has described the metheds used to 
design individual blade sections. When iocated as desired 
relative to the blade stacking axis, the section coordinates are 
the coordinates of the streamsurface blade section. A series 
of sections on ali streamsurfaces specifies the envelope of the 
blade, but the surface coordinates are not in a form convenient 
for manufacturing purposes. The calculation procedure uses a 
spline-curve to interpolate (or extrapolate) the coordinates of | 
the blade surfaces for manufacturing purposes. 


3. SECTION PROPERTIES 


The staeking axis of the blade is passed through each 
streamsurface section either at one of the edges ov at a point 
specified relative to the centroid of the section. Because the 
streamsurface sections are in general non-planar, the centroids 
of the manufacturing sections will not generally lie precisely 
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om the stacking axis when the streamsurface sections are stacked 
on their centroids. By determining the locations of the 
centroids of the manufacturing sections so obtai' ed, it is 
possible to estimate the offsets that must be applied wher 
restacking the streamsurface sections to locate the manufacturing 
centroids as desired relative to the stacking axis. 


To assist further in the mechanical analysis cf the blade, 
the areas, second moments of area, principal axes, and principal 
second moments of area for both the streamsurface and manu- 
facturing sections are also determired in the calculation 
procedure. 


4, BLADE CHARACTERISTICS 


& calculation of the volume enclosed by the blade between 
the innermost and outermost streamsurfaces is made. In addition, 
quantities which describe the blade on cylindrical surfaces and 
which may be required in an aerodynamic analysis of the blade 
are computed as an option. The calculations are presented here 
because a typographical error undetected during editing of 
Reference 1 has impaired their usefulness in that Reference. 


First, the angular position of the camber line with respect 
to the stack axis at a streamsurface-computing station 
intersection is specified in terms of %, defined in Figure l. 


The physical passage blockage (B) due to the presence of 
the blades is determined as a percentage of the passage 
circumference in terms of the number of blades in the blade 
row and T, the angie subtended on the cylindrical surface by 
each blade: 


B= Nt 


a (25) 


The blade lean angle, © , with respect to the radial 
direction at a given point is obtained from the slope of a 
spline~curve fit through the y-z Cartesian coordinates of the 
streamsurface section camber lines at the particular axial 
location. 


Thus 


€= m - Arctan (2x) (26) 
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Two other quantities are needed to produce the proper 
méan-camber line angle cn the cylindrical surface: The local 
computing station inclination, yw, obtained from the specified 
Station description; and the local Streamsurface inclination, y, 
obtained from the specified x-r streamsurface description. 
Together with the camber line angle on the streamsurface, 
these quantities are employed in the following equation to 
calculate the proper cylindrical-surface section angle, ag: 
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tan y tan € + tan a,/cos y 


1 - tan uw tan vy 


tan aZ = 


(27) 
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‘SECTION V = 


USE OF THE PROGRAN 


Basic information required by the user to run the ARL 
computer program which incorporates the calculation procedure 
described in this report is given in this section. The various 
input data items are defined first, and the input data rormat 
is then specified. A description of the output data that may 
be expected is given. (Implementation of the program on a 
computing system is not discussed here. but in the section 
entitled "Computer Program Details". ) 


1. DEFINITION OF INPUT DATA ITEMS 


TITLE An alphanumeric title of 72 characters that may be 
used to identify a run. 


NLINES The number of streamsurfaces which are defined and 
on which blade sections will be designed. Must 
satisfy 2<NLINES <15. 


NSTNS The number of computing stations at which the 
streamsurface radii are specified. Must satisfy 
3S NSTNS < 10. 


NZ The number of constant-z planes on which manu- | 
facturing (Cartesian) coordinates for the blade 
are required. Must satisfy 3<NZ2<15. 


NSPEC The number of radially-disposed points at which 
the parameters of the blade sections are specified. 
Must satisfy l1< NSPEC S15. 
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ISEGPT The number of points to be used to define each an : ! 

1. . < Sete : 

segment of the camber line 2 < ISEGPT < Integer (rapg tae , 

NBLADE The number of blades in the blade row. : 
ISTAK If ISTAK = 0, the blade wiil be stacked at the 


leading edge. 


If ISTAK = 1, the blade will be stacked at the 
trailing edge. 


If ISTAK = 2, the blade will be stacked at, or 
offset from, the section centroid. 


LY 


IPUNCH 


LIFPLOT 


IPRINT 


ZINNER, 
ZOUTER 


SCALE 


If TPUNCH = 0, the quantities necessary for aero~ . 
dynamic analysis of the resulting blade are not ; 
preduced on punched cards. 


If IPUNCH = 1, these quantities are produced on 
punched cards. 


Where CALCOMP software is incorporated into the 
computing system, IFPLCT specifies the creation of 
precision plots. (Further information regarding the 
requirements for this are given in the section 
entitled "Computer Program Details.")- 


If IFPLOT = 0, no plots will be produced. 


If IFPLOT = 1, a plot of the streamsurface sections 
Will be preduced. 411 NLINES sectitns are shown 
superimposed. The origin for each section plot is 
offset from the centroid of the section by distances 
specified by DELX and DFLY. If IFPLOT = 2, a plot 
of the manufacturing sections will be produced. The 
origin is the blade stacking axis, and all NZ 
sections are shown superimposed. 


If IFPLOT = 3, both of the plots described for 
IFPLOT = 1 and 2 will be produced. 


If IFPLOT = 4, individual plots of each of the 

manufacturing sections will be produced. The axes | 
are rotated clockwise by the section stagger angle 3 
for each plot. | 


The input data is always listed by the program. 
Details of the streamsurface and manufacturing 
sections are printed as prescribed by IPRINT. 


If IPRINT = 0, details of streamsurface and 
manufacturing sections are printed. 


If IPRINT = 1, details of streamsurface sections 
are printed. 


If IPRINT = 2, details of the manufacturing 
sections are prixced. 


The NZ manufacturing sections are equispaced between 


Zz equals ZINNER and ZOUTER. ; 


When precision plots are produced, SCALE is the scale 
factor employed. 
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NEM SA re 


STACKZ 


PLTSZE 
ERLE 
IRT&S 
NRADEV 
NINC 


NSIGN 


IFCA 


IPASS 
AKSHPE 
SOLTOL 
NPTS 
RADEV 


OM 


DEVCRYV 


hae ee —~_~ -- Fe ~_Ss 


The axial coordinate of the stacking axis for the 
blade, relative to the same origin as used for the 
station locations, XSTA. 


The size (inches) of the plotter to be used in the 
creation of precision plots. 


The number of the computing station designated as 
the blade leading edge. 


The number of the computing station designated as 
the blade trailing edge. 


The number of radii at which the non-dimensional 
deviation distribution is specified. 1S NRADEVS 5. 


The number of points which describe the incidence 
angle distribution. 1<NINC S15. 


An integer which specifies the sign convention of 
the particular blade. Conventionally positive 
rotors and stators have NSIGN values of -1 and +1, 
respectively. 


If IFCA = 1, the factor m in the deviation angle 
rule is that of the NACA-~-65-series mean line 
(Figure 195, Reference 2). 


If IFCA = 2, the factor m in the deviation angle 
rule is that of the circular-are mean line. 


The number of initial values of S/R, which are to 
be used in the procedure to find the optimal 
camber line. 20S IPASS < 50, 


the shape factor (Kg .) in the deviation equation. 


The solidity tolerance used in the iterative 
procedure to produce a consistent camber line, 


The number of points defining a particular chord- 
wise deviation distribution. lSNPTS<10. 


Radius at which a particular deviation distribution 
applies. 


An array of NPTS meridional chord fractions which, 
together with DEVCRV, specify a particular 
deviation distribution, 


An array of NPTS normalized deviation fractions which, 


together with SM, specify a particular deviation 
distribution. 
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RINC 


XINC 


DELDEV 


KPTS 


IFANGS 


XSTA 


AITRANG 


An array of NINC radii which, together with XINC, 
specify the incidence distribution at the leading 
edge. 


An array of NINC incidence angles which, together 

with RINC, specify the incidence distribution. 

Input positive for cenventionally positive rotors 

and stators (see NSIGN). 


An array of NINC angles which, together with RINC, 
specify the distribution of the “arbitrary extra 
deviation" term in the deviation determination. 
Input positive for conventionally positive rotors 
and stators (see NST). 


The number of points provided to specify the shape 
of a computing station. 


If KPTS = 1, the computing station is upright and 
linear. 


If KPTS = 2, the computing station is linear and 
either upright or inclined. 


If KPTS > 2, a spiine curve is fitted through the 
points provided to specify the shape of the station. 


If ITVANGS = 0, the caiculations of the quantities 
required for aerodynamic analysis will be omitted 
at a particular computing station. 


If IFANGS = 1, these calculations will be performed 
at that station. 


An array of KPTS axial coordinates (relative to an 
arbitrary origin) which, together with RSTA, specify 
the shape of a particular computing station. 


An array of KPTS radii which, together with XSTA, 
specify the shape of a particular computing station. 


The streamsurface radii at NLINES locations at each 
of the NSTNS stations. 


The relative air angles at NLINES locations at each 
of the NSTNS stations. 
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ZR The variation of properties of the streamsurface 
blade sections is specified as a functicn of 
streamsurface number. The various quantities are 
then interpclated (or extrapolated) at each strean- 
surface. The streamsurfaces ere numbered 
consecutively from the innermost ontward, starting 
with 1.0. ZR must increase monotonically, there 
being NSPEC vaines in all. 


YA The fraction of meridional chord used as the leading 
edge in the calculation of the section chord for 
the solidity iteration on 2 particular streamsurface. 
If YA = 0., the true chord length is calculated. 
O.0SYA<1.0. 


Y3 The increment in S/Rp which, with SDIVYR and IPASS, 
establishes the initial range of S/R, which is 
inspected in the determination of the optimal 
camber line on a particular streansurface. May be 
positive or negative. 


ye If YC = 0., the radius of curvature at the leading 
edge of each camber line will be considered in the 
procedure to identify the camber line which 
maximizes the minimum radius of curvature. 


If YC = 1.0, the vadius of curvature at the leading 
edge wiil not be considered in this procedure. 


YE The maximum number of inflection points expected on 
a particular camber line. If the calculated 
minimum number is greater than YE, an informational 
diagnostic is printed. 


RLE The ratio of section leading edge radius to chord. 

TC The ratio of section maximum thickness to  hord. 

TE The ratio of section trailing edge half~—thickness 
to chord. 

Z The location of the section maximum thickness, as ,a 


fraction of camber line length. 


SDIVR The initial value of S/R,. 
DELX, The stacking axis passes through the streamsurface 
DELY blade sections, offset from the centroid, leading, 


or trailing edge by DELX and DELY in the x and y 
directions, respectively. 
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2. <INPUT DATA FORMAT 

Data input is by punched card, and three formats are used. 
The first card oniy is alphanumeric, using the first 72 coiumns 
of the card. Integers are placed in three-cclumn fields, which 
start with Column 1. No decimal points are used, and the 
integer should be right-justified. Real numbers are placed in 
12-column fields, which also start with Column 1. Decimal 
points should be included, and the numbers may be placed 
anywhere in the field. 

In the following chart, one line corresponds to one card. 
TITLE 
NLINES NSTNS NZ NSPEC ISEGPT NBLADE ISTAK IPUNCH IFPLOT IPRINT 
ZIHNER ZOUTER SCALE STACKX PLTSZE 
IRLE IRTE NRADEV NINC NSIGN IFCA IPASS 
XKSHPE SOLTOL 
NPTS 
RADEV repeated NRADEV times | 
SM DEVCRV } repeated NPTS times 
RINC XINC DELDEV } repeated NINC times 
KPTS IFANGS 
XSTA RSTA } repeated KPTS times repeated NSTNS times 
R ALRANG } repeated NLINES times 
ZR YA YB YC YE RLE } ! 

f repeated NSPEC times 

TC TE %2 SDIVR DELX DELY 


Listing of a sample input data deck is included under 
"Example of Use of the Program." 


3. OUTPUT DATA 


Printed output from the program may be considered to 
consist of four sections: a printout of the input data, details 
of the camber line and blade section on each streamsurface, a 
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listing of quantities required for aerodynamic analysis, and 
Getails of the manufecturing sections determined on the 
constant-Z planes. These are briefly described below. 


The input data printout includes ail quantitits read ia, 
and is self-explanatory. 


Details of the streamsurface blade sections are printed if 
ITPRINT = Q or 1. Listed first are the results of the 
investigations of the S/Rp parameter. The initial table presents 
the results of the first iteration which is used to idenzity 
the range of S/R, in which the minimum number of inflection 
points occur on the camber line. This range is in tura 
investigated with finer increments of S/Rp to determine the 
maximum value of the minimum radius of curvature. Then follow 
the details of the optimal camber line which has been identified 
by a third investigation of S/Rp with still finer parameter 
increments. These details include the deviation and solidity 
which have been calculated for this optimal S/Ro, ani a 
description of the camber line in terms of coordina*es, first 
and second derivatives, and the radii of curvature. Listed next 
are the parameters defining the blade section, some of which 
are computed and some of which are interpolated ai the stream- 
surface from the tables read in. Then follow details of the 
blade section in "normalized" form. The blade section geometry 
is given for the particular section, except that the meridional 
projection of the chord is unity. For this section of the 
output, the coordinate origin is the blade lead‘ng edge. The 
following quantities are given: blade chord, stagger angle, 
camber angle, section area, location of centroid of the section, 
second moments of area of the section about the centroid, 
orientation of the principal axes, and the principal second 
moments of area of the section about the cenlroid. Then are 
listed the coordinates of the camber line, the camber line angle, 
the section thickness, and the coordinates cf the blade surfaces. 
A line-printer plot of the normalized section follows. The 
scales for the plot are arranged so that the section just fills 
the page. Thus, the scales will generally differ from one plot 
to another. "Dimensional" details of the blade section are 
given next. The normalized data given previously is scaled to 
give the proper biade section. For this section of the output, 
the coordinates are with respect to the blade stacking axis. 

The following quantities are given: blade chord, radius and 
location of center of the leading edge, section area, the 

second moments of area of the section about the centroid, and 

the principal second moments of area of the section about the 
centroid. The coordinates of points on the blade surfaces are 
then listed, followed by the coordinates of 31 points distributed 
at six degree intervals around the leading edge. Finally, the 
coordinates of the blade surfaces and points around the leading 
edge are shown in Cartesian form. 
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The quantities required for aerodynamic analysis are 
printed at ail computing stations specified by the IFANGS 
parameter. The radius, section angle, blade lean angie, blade 
blockage, and relative angular location of the camber line are 
printed at each streamsurface intersection with the particular 
conputing station. 


Details of the manufacturing sections are printed if 
ITPRINT = 0 or 2. At each value of z specified by ZINNER, 
ZOUTER and NZ, section properties and ccordinates are given. 
The origin for the coordinates is the blade stacking axis. The 
following quantities are given: section area, the location of 
the centroid of the section, the second moments of area of 
the section about the centroid, the principal second moments of 
area of the section about the centroid, the orientation of the 
principal axes, and the section torsional constant. Then the 
coordinates of points on the blade section surfaces are listed, 
followed by 31 points around the leading edge. 


Precision piots are produced if IFPLOT = 1, 2, 3 or 4 as 
described under the definition of £FPLOT given previously. 


If TPUNCH = 1, the program punches the quantities required 
for aerodynamic analysis, together with identifying indices 
denoting station number and streamsurface number, on cards in 
the following format: 5 fields each of 12 locations for the 
quantities themselves, followed by 2 fields each of 3 locations 
for the indices. 


ey ee 
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SECTION Vi 
EXAHPLE OF USE OF THE PROGRAM 


This section shows the use of the program to generate a 
compressor rotor with an inlet hub-to-tip ratio of 0.31, a rather 
steep hub ramp angle (32.5°), and a constant outer radius. The 
blade is defined by six computing stations, one at either edge 
and four internal. This results in a camber line composed. of 
five segments. Six streamlines have been used to define the 
flow and hence the blade by means of the streamsurface blade 
sections. The computing-stations and streamsurfaces which 
aefine the blade are depicted in a stack-axis projection in 
Figure 2. 


1. INPUT DATA 


The input data deck used for this example is listed below. 
Some ~oints of interest are noted in the order in which they 
occur in the input. 


As mentioned above, six streamsurface bladé sections are 
used to define the blade. The streamsurface radii are specified 
at eight computing stations, the first ani last of which are 
outside of the rotor. This ensures that the boundary conditions 
imposed on the spline-curve (zero curvature at the endpoints) 
have little influence on the shape of the curve representing 
the streamsurface within the blade. The useful relative air 
angles are specified at the six computing stations defining the 
blade. The computing stations within the blade are curved in 
an attempt to make the meridional projections of the camber line 
segments approximately equal. The parameters which define the 
streamsurface blade sections are given at six locations; that 
is, at every streamsurface. Twelve points are used to define 
each segment, which results in a camber line defined by a total 
of fifty-six points. (This number serves the purposes of this 
example well, but it would be advantageous to use more points 
if the precision-plot output is to be incorporated directly in 
the manufacturing proc:dure.) There are twenty blades in the 
particular blade row, stacked approximately on their stream- 
surface centroids. No punched output is requested. All 
optional sections of the printed output are to be printed, and 
superimposed section plots are to be produced on a plotter 
(with an 1ll-inch useful range) at two and one-half times fuil 
size. 


Stations 2 and 7 are the leading and trailing edges of the 
blade, respectively. The deviation distribution is the same 
at all radii, since it is specified only at one radius. The 
incidence angles and the arbitrary extra deviation are specified 
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at six points for this votor blade. The factor m in the deviation 
caiculation is to be that of a circuvlar-arc camber line. Thirty 
camber lines will initially be investigated in the erfort to 
identify the optimal camber line, representing a broad range of 
the S/R, paraneter. The shape factor in the deviation 

calculation is unity, and the tolerance on solidity in the 
iteration for each camber line is 0.005. 


The streamsurface radii and relative air angles have teen > 
leterrined from an aerodynamic analysis of the flow through the 
bladé row. The leading edge of the blade has been established 
as the point from which the chord length, and thus the solidity, 
will be computed. 


The increment in the S/R, parameter has been initially set 
equal to + 0.04 in the investigation to find the optimal camber 
line. This increment will autometically be reduced twicz 
subsequently as the most favorable range of S/Rp is nore 
closely scrutinized. It is anticipated that no inflect on 
points will be required near the hub, while a single “nflection 
point wilJ. probabiy be required in the mid and tip re ions of 
the blade. The blade thickness distribution is detertined by 
aerodynamic and mechanical factors. The leading edge radius 
and trailing edge half-thickness are set tc approximately 0.005 
inch, probably a practical minimum. Maximum thickness/ chord 
ratios vary from 6 percent at the hub to 2.5 percent at the 
casing. The maximum thickness is placed in the rearward portion 
of the camber line, which helps to maintain a small leading 
edge wedge angle. 
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2. QUTPUT DATA 


The input data specified ali optional output data, and a 
Sample of each segment of the output is presented in this 
section. A brief description of some espects of the output is 
presented below. . 


Shown first is the printout of the input data. This is 
Followed by details of the First streamsurface blade section. 
The tables presenting the results of the S/R, investigations are 
followed by a description of the camber line which has been 
selected. Then appear the details of, and a line-~printer plot 
of, the normalized blade section, followed by the specifications 
of the section scaled ton the proper dimensions. Printed first 
are the 56 points specified for each blade surface. Next 
appear the 31 points caescribing the leading edge radius. The 
Final data for the strearsurface section are the equivalent 
Cartesian coordinates for these same points. 


The format is repeated for each of the remaining stream- 
surface blade sections, but these results are not reproduced 
here. Subsequent to the streamsurface data are printed the 
quantities required for aerodynamic analysis at those stations 
where requested by the IFANGS parameter. 


Details for the seven manufacturing sections defining the 
blade follow. Reproduced below is the output relating to the 
first (innermost) section. Properties of the section are 
followed Ly coordinates of 56 points on each surface and 31 
points around the leading edge. It will be noted that the 
section centroid is not calculated to be exactly on the stacking 
axis. If it were desired to have the centroids of the manu- 
facturing sections lie more nearly precisely on the stacking 
axis, the program would be rerun with either DELX and DELY 
offsets specified so that they would counteract the mislocation 
of the centroids previously determined, or with a slightly 
shifted location of the stack axis. However, the user must bear 
in mind that the final step in the design procedure is to obtain 
Satisfactory coincidence of the blade with the stations defining 
its leading and trailing edges in the meridional plane. It is 
possible that a blade which had all manufacturing sections 
stacked precisely on their centroids at a particular stacking- 
axis location would have quite an undesirable meridional profile. 
There are tradeoffs, then, between the stacking preciseness and 
the coincidence of the calculated blade's edges with its assumed 
profile on the one hand, and the number of iterations it would 
require to find the truly optimai stack location and centroid 
offsets, on the other. 


The input data specified superimposed precision plots of 


both the streamsurfdce sections and the manufacturing sections. 
These plots are reproduced (at reduced size) as Figures 3 and 4, 
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respectively. It is of interest to refer also to Figure 2. The 
innermost manufacturing plane is well below the lowest point on 
the hub streamsurface section in the plane of the stacking axis. 
The Cartesian coordinates of this streamsurface section show 
that the lowest pcint on the section (the 14th point on the 
leading edge) is at Z = 2.5722h, The streamsurface radius at 
this point is 2.6514, and the innermost manufacturing plane is 
at Z= 2.5. Thus, at the leading edge, the extrapolation 
requirce to define the manufacturing section is somewhat smaller 
than might first appear. t the trailing edge, extrapolation 

is required for the first three manufacturing sections. The 
streamsurface radius at the casing and the Z-coordinate for 

the outermost mauufacturing plane both equal 8.53; however, the 
Z-coordinates of the blade section are all actually below 8,5. 
Thus, the outermost section too is defined completely by extrap- 
olation. Of course, portions of the blade that are defined by 
extrapolation do not appear on the final blade, but facilitate 
manufacture. 
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SECTION Vil 
COMPUTER PROGRAM DETAILS 


1. FMPLEMENTATION OF THE COMPUTER PROGRAM 


The program is written in FORTRAN IV and was developed on 
a CbdC 6000 Series System, incorporating the version 3.3 SCOPE 
Operating System. When loaded into core, the program (and 
resident system) occupies about 32K of storage, so that the 
program will probably not be usable without modification on a 
relatively small computer. Apart from this limitation, the 
program should be compatible with the majority of modern 
computing systems. The program consists of .a main program, and 
Subroutines BQ, CQ, Dl, EQ, and FQ. The six decks have been 
given the identifiers A, B, C, D, E, and F, respectively. 
Listings of the decks are shown beiow, and the deck set-up 
required for the CDC 6000 Series System is also presented. 


The program uses three numerical system units for its 
input and output routines. Input is drawn from the card unit 
via READ statements referring to Unit LOG 15 output is sent to 
the line printer by WRITE statements referring to Unit LOG 2; 
and punched output is produced via WRITE statements referenced 
to Unit LOG 3. Units LOG i, LOG 2, and LOG 3 are set equal to 
5, 6, and 7, respectively, on cards A130-140 in the FORTRAN 
programming. On the CDC 6000 Series System, the "PROGRAM" card 
must also establish the input and output linkages. On other 
computing systems, the "PRGGRAM" card may not be required, and 
the input-output files may be established via control cards. 


The program as presented herein utilizes an on-line 
precision plotting capability available at the program 
development site. Calls to four subroutines not included in 
the deck are included in the program. These calls are executed 
only if precision plots are specified in the input data, and 
the entry points expressed are part cf standard CALCOMP software 
normally supplied to users of CALCOMP precision plotting 
equipment. The various call statements used in the program are 
explained below, to facilitate modification should the need arise. 


CALL PLOT (XPLOT, YPLOT, N) 


The majority of plotting is done using this form of call. The 
parameters XPLOT and YPLOT are the "x" and "y" coordinates (in 
inches) on the paper to which the pen is being directed. The 
parameter N indicates pen up or down, | 1 | = 3 or |N] = 2, 
respectively, and will cause XPLOT or YPLOT toa be assigned 

as the origin for further coordinates if N is negative. 
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CALL SYMBOL (X, Y, H, TEXT, THETA, 2B) 


This call is used to title the plots. The parameters X and Y 
are the coordixates (in inches) of the Lower left hang corner 
of the first character, H is the character height (in inches), 
TEXT is the character to be printed, THETA is the angle of the 
lettering with respect to the "x" axis and N is the total 
number of characters to be printed. 


CALL NUMBER (X, Y, H, F, THETA, N) 
This cal] causes tne printing of the number F. The parameters 
X, ¥, H, and THETA are used as for CALL SYMBOL. The parameter 


N indicates the number of digits following the decimal point 
if positive, or truncation to an integer if equal to -1l. 


CALL PLOTE 
This call terminates the tape. 

In the event the program is used on a computing system 
which does not include CALCOMP software, and the operating 
system wili not execute a program with unsatisfied external 
references, dummy entry points may be supplied by adding to 
the deck a subroutine such as the following: 

SUBROUTINE PLOT 

A= A ; 
ENTRY SYMBOL 

ENTRY NUMBER 

ENTRY PLOTE 

RETURN 


END 


z2. DECK SETUP FOR CDC 6000 SERIES SYSTEM 


ihe deck setup required to run the program on a CDC 6000 
Series System incorporating the SCOPE 3.3 Operating System is 
shown below. Production runs of the program would usually 
employ relocatable binary forms of the routines produced from 
the source decks to avoid having to compile the FORTRAN for 
each run and the waste of associated computer resources. 
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JOB Identification, etc. 


FTN. 

LGO. 

7/8/23 End of Record 
SOURCE DECK A 
SOURCE DECK B 
SOURCE DECK C 
SOURCE DECK D 
SOURCE DECK E 
SOURCE DECK F 

7/8/83 End of Record 
DATA DECK 


6/7/8/9 End of Job 


3. FORTRAN PROGRAM LISTING 


A listing of the FORTRAN program appears on the following 
pages with each subroutine started on a new page. 
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4. PROGRAM LOGIC 


The calculation procedure which has been described in 
this report is performed primarily in the main program and 
Subroutine 8Q. The calculations regarding the construction of 
the camber line, the stacking of the streamsurface blade 
sections, and the determination of manufacturing sections are 
performed in the main program. Subroutine BQ applies the 
thickness distribution to the camber line and determines 
quantities necessary to obtain the Cartesian coordinates of 
the section from the streamsurface coordinates. Subroutine D1 
is the curve-fitting routine, and CQ is used to determine 
slopes of varicus spiine curves at particular points. EQ 
produces the line-printer section plot in the printed portion 
of the output, and FQ produces axes on the section plots for 
IFPLOT = 1. 2, or 3. 


& description of the calculation procedure employed in 
the main program and in Subroutine BQ is described below. 
Each step is keyed to its location in the program by the 
parenthetical deck serialization. 


1. The input data is read and printed. (A155-A680) 


2. If precision plotting is specified, the plot is initialized, 
and axes produced if IFPLCOT = 1 or 3. (A695-A710) 


3. <A loop which creates a section on each streamsurface is 
commenced. (A715) 


4. The axial locations of the intersections of a particular 
streamsurface with the computing stations describing the blade 
are determined. The meridional streamsurface length is 
obtained as described in Equation (1). (A720-A780) 


5. The parameters relating to the streamsurface blade section 
are interpolated (or extrapolated) from the input tables. If 
NSPEC = 1, they are taken to be radially uniform. If NSPEC = 2, 
linear interpolation is used; if NSPEC = 3, spline-curve 
interpolation is employed. (A785-A840) 


6. The loop to determine the optimal camber line is initiated. 
(A855) 


7. The first estimate of true chord is calculated per 
quation (2), and the solidity as in Equation (3). (A870~A910) 


8, The incidence angle and extra deviation applicable to the 

particular section are obtained by interpolation of the radius 
at the leading edge of the streamsurface in the input tables. 

(A915-A920} 
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S. The section angle at the leading edge is determined as in 
Equation (6). (A925) 


10. The quantities reaquirei for the deviation angle caiculatioa 
are obtained by interpolation from various figures of 
Reference 2. (A935-A960) 


ll. The deviation angle is calculated using Equation (5). 
(A995-A1000) 


12. The section angie at the trailing edge is calculated using 
Equation (7), and at internal points using Equation (8), with 
fractions of deviation obtained by radial interpolation from 
the input distributions based on the streamsurface radius at 
the trailing edge of the blade. (A1020-A1085,) 


13. A camber line is constructed of cubic segments following 
the analysis of Equations (9) - (15) for the initial value of 
S/Ry- The number of inflection points is determined. (A1090- 
A1350)} 


14. The iteration on solidity is performed until the tolerance 
is within the prescribed limit. (A1355-A1395) 


15. Steps 13 and 14 are repeated for IPASS-1 values of the 
S/Ro parameter. (A140G-A1450) 


16. The range of S/Ro with the minimum number of inflection 
points is established.: (A1480-A1600) 


17. Steps 13, 14 and 15 are repeated for finer increments of 
the S/R, parameter in the range determined in 16. The maxinum 
value of the minimum radius of curvature in this range is 
determined. {(A1645-A1690) 


18. Stili finer increments of S/R, on either side of the S/Ro 
value which produced the maximum value of the minimum radius of 
curvature in 17 are examined to find the optimal S/Ro value. 
(A1730-A1745) 


19. The details of the optimal camber line are printed if 
IPRINT = 0 or 1. (A1755-A1810) 


20. The normalized chord length of the optimal camber line is 
computed. (A1820) 


21. The total streamsurface length is calculated. (A1830-A1890) 


22. If IPRINT = 0 or 1, the parameters defining the section 
are printed. (B65-B115) 
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23. The coefficients of the twe thickness equations are 
computed. (B175-5220) 


24%, At each point of the camber line, the corresponding 
coordinates on each bliade surface are obtained from the 
coordinates and slope of the camber line, and the appropriate 
thickness, scaled for an overall section meridicnal chord of 
unity. (B225-B310) 


25. The section 2rea and centroid location are determined. 
(B315-B365) 


26. The camber line is redefined in terms of 100 points to 
assure sufficient points fer accurate linear interpolation in 
the determination of 9 (Figure 1), needed for the eventual 
Cartesian coordinate determination. (B380-B440) 


27. The various streamsurface section properties are 
determined. (B545-B560) 


28. If IPRINT = 0 or 1, details of the normalized blade 
section and a line-printer plot are produced. (B575-B750) 


29. Sectional properties are scaled to produce the "dimensional” 
results. (B755~B830) 


30. If IFRINT = 0 or 1, the section information is printed. 
(B840-B1095) 


31. The coordinates of 31 points describing the leading edge 
are determined. (B1070~B1080) 


32, The origin of the coordinate system is shifted to the 
stack axis, and the relative ¢@ values for the blade surfaces 
are determined. (B1120-B1595) 


33. If IFPLOT = 1 or 3, the streamsurface section plot is 
produced. (A1950-A2040) 


34. If the calculations for aerodynamic analysis are required, 
various items related to the camber line are stored. (A2045- 
A2095) 


35. The Cartesian coordinates for each point on the section 
surface are computed, and printed if IPRINT = 0 or 1. (A2100- 
A2435) 


36. The loop initiated in 3 is repeated for each streamsurface. 
(A2440) 


37. Unless IPRINT = 1, the volume of the blade is computed and 
printed. (A2445-A2555) 
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38. If specified, the calculations for zeredynanic analysis 
are performed and printed, and punched if IPUNCH = 3. (A2560- 
A2790) 


39. if IFPLCT = 2 or 3, the exes are drawn and titled for the 
superimposed plot of the manufacturing sections. (A2806) 


G9. If no output relating to manufacturing sections is specified 
by either IFPLOT or IPRINT, the remainder of the progran is 
bypassed. Alternatively, if printed detaiis of the rmanufactur- 
ing ne are specified, a heading is printed. (A2895- 

A2335 


41, “he location of each of the manufacturing planes is 
deternined. (A2840-A2865) 


42. The (Cartesian) coordinates of each point on the biade 
surface are obtained by spline-curve interpolation at each of 
the manufacturing sections. (A2870-A2945) 


“3. A loop that is performed for each manufacturing section is 
initiated. This loop contains the determination of section 
properties and the output of results for the section. If 

ITPRINT = 0 or 2, section properties and coordinates are printed. 
(A2950-A3495) 


ua. If IFPLOT = 2 or 3, a plot of the manufacturing sections 
is produced. (A3510-A3605) 


G5. If IFPLOT = 4, an individual plot of the manufacturing 
sections is made. The axes are rotated clockwise until the 
chord line is horizontal. The angle cf rotation is indicated 
as the stagger angle. (A3610-A3845) 


46. The’ loop initiated in Step 43 for each manufacturing 
section is terminated. (43850) 


u7. If precision plots have been made, the plotting is 
terminated. (A3855) 
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Example Blade Design: 
Streamsurface Sections 
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